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We present experiments on the magnetic held-dependent thermal transport in the spin-1/2 ladder 
system (C5Hi2N)2CuBr4. The thermal conductivity k.{B) is only weakly affected by the field- 
induced transitions between the gapless Luttinger-liquid state realized for Bd < B < Bc2 and the 
gapped states, suggesting the absence of a direct contribution of the spin excitations to the heat 
transport. We observe, however, that the thermal conductivity is strongly suppressed by the mag- 
netic field deeply within the Luttinger-liquid state. These surprising observations are discussed in 
terms of localization of spinous within finite ladder segments and spinon-phonon umklapp scattering 
of the predominantly phononic heat transport. 



Studies of the heat transport in ID spin systems 
are of strong current interest. [TJ [51 |3] From the the- 
oretical side there is consensus that the intrinsic spin- 
mediated heat transport of integrable spin models is bal- 
listic, while the situation in nonintegrable spin models 
is less clear. Experimentally, such studies were stimu- 
lated by the observation of a strong anisotropy of the 
thermal conductivity k(T) in the spin-1/2 ladder com- 
pound (Sr,Ca,La)i4Cu2404i, which has been explained 
by a large spin contribution Ks along the ladder direc- 
tion [H m H]. In order to relate model calculations to 
experimental data, theory has to incorporate the cou- 
pling between spin excitations and the underlying lattice, 
while experimentally it is necessary to separate Ks from 
the measured total k. Here, studies of the magnetic- 
field dependent n{B, T) can provide much more informa- 
tion than just the zero-field k(T), since strong enough 
magnetic fields change the spin excitation spectra and 
cause transitions between different quantum phases. Yet, 
the above-mentioned (Sr,Ca,La)i4Cu2404i is unsuitable 
for such studies due to the strong exchange interaction 
(~ 2000 K) in the ladders that requires magnetic fields 
far above typical laboratory magnets. 

Piperidinium copper bromide (C5Hi2N)2CuBr4 [7J is 
one of the best spin-1/2 ladder compounds with weak in- 
traladder exchange. It has a monoclinic structure (space 
group P2i/c) with the ladders running along the a axis. 
The rung exchange Jj^ ~ 13 K is about four times larger 
than the leg exchange Jy ~ 3.6 K [iillTniinilllllllllll- 
The zero-field spin excitation gap A ~ 9.5 K decreases in 
an external magnetic field, so that between Bd ~ 7.1 T 
and the full saturation field Bc2 — 14.5 T (for B \\ a 
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where g = 2.06 |7]) the spin excitations are gapless and 
the system is in the Luttinger liquid (LL) state. The 
LL state extends down to the temperature of a 3D or- 
dering Tn < 110 mK suggesting an interladder coupling 
J' -27mK<C J_L, J|| [iniinilTS]. 

Here, we present measurements of the thermal con- 
ductivity of (C5Hi2N)2CuBr4 between 0.3 and 10 K in 
magnetic fields up to 17 T. Surprisingly, no sign of spin- 
mediated heat transport is observed, which is interpreted 
by spinon localization. The heat transport is dominated 
by phonons being strongly scattered by spin excitations 
in the vicinity of the commensurate wavevector. 

Two samples of (C5Hi2N)2CuBr4 of approximate di- 
mensions 5 X 1.7 X 0.8 mm'^ were cut from two crystals of 
the same batch with the longest dimension either along 
the a axis (along the ladders) or along the c* axis. The 
thermal conductivity was measured using the standard 
uniaxial heat-flow method where the temperature dif- 
ference was produced by a heater attached to one end 
of the sample and monitored using a matched pair of 
Ru02 thermometers. The heat flow was directed along 
the longest direction of each sample and the magnetic 
field was parallel to a. 

The zero-field thermal conductivities k" , k-^ for both 
directions of the heat fiow are shown in Fig.[l](a). Their 
behavior is typical for a phononic heat transport, ap- 
proaching K cx T'^ at lowest temperature. The field de- 
pendencies of the thermal conductivity normalized to its 
zero-field value at several constant temperatures are pre- 
sented in Fig. [l] (b). The most pronounced feature of 
the k{B) curves are two minima located in the middle 
of the field interval between Bd and Bc2 symmetrically 
with respect to Bmiddio = (Bd + Bc2)/2 ~ 10.8 T. With 
decreasing temperature, the two minima become deeper 
and move closer to i?middic- It is noteworthy that the 
critical fields Bd and i?c2 are not marked by a distinct 
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FIG. 1: (color online) (a ) Zero-field thermal conductivity 
of (C5Hi2N)2CuBr4 parallel and perpendicular to the ladder 
direction as a function of temperature, (b) Relative change of 
the thermal conductivity parallel to the ladders as a function 
of magnetic field at constant temperatures. Inset: Data for 
T = 0.77 K parallel and perpendicular to the ladders. 



feature of k{B). As displayed in the inset of Fig. [T] (b), 
kII(B) and K-^{B) show the same behavior. 

The total heat conductivity of a magnetic insulator 
comprises the phonon and the spin contributions 

K = Ks + Kph, (1) 

where each contribution is given, as a first approxima- 
tion, by a product of the specific heat Ci, the average ve- 
locity Vi, and the mean free path £i of phonons (i ~ ph) 
and spin excitations (i = s). In addition, a term in 
Eq. ([T]) due to spin-phonon drag is expected [TS]. In 
quasi-lD systems, both Ks and spin-phonon drag are only 
essential along the chain (ladder) direction. From the al- 



most absent anisotropy of k}^{B) and k-^{B) we conclude 
that Ks and the drag contribution are small. 

The spin contribution k,s{T,B) can be analyzed using 
a mapping of the spin-1/2 ladder onto the effective spin- 
1/2 XXZ chain. Within this mapping, the ground state 
singlet and the low-energy component of the triplet form 
an effective spin S = 1/2 on each rung of a ladder [THJ 
[T7| [T8] . The effective Hamiltonian is 

N 

Hxxz = jY.(si~Sl+'+~Sl~Sl+'+5~Sl~Sl+') 

i 

N 

- g^iBBJ2Sl, (2) 

i 

where S is the anisotropy parameter. For a spin-1/2 lad- 
der in the strong-coupling limit {J±/J\\ ^ 1), 6 — 1/2, 
J = J||, and S = S - {Jj_ + J||/2)/g/ii3. Thus, B = 
^middle corrcspouds to the zero effective field B — 0. The 
model hamiltonian ^ well describes many observed low- 
temperature features of (C5Hi2N)2CuBr4, as e. g. the 
magnetization [H] , thermal expansion [9] , NMR [TD] , and 
specific heat |12j measurements. Some discrepancies are 
mainly caused by the fact that the condition J±/J\\ » 1 
is only approximately satisfied. 

The heat transport in the spin-1/2 XXZ chain can be 
well described for T < J/k^ within a relaxation time ap- 
proximation combined with a mean-field theory (MFT) 
approach via the Jordan- Wigner (JW) transformation. 
Eq. (|2| is mapped onto a system of interacting spinless 
fermions |191 120j , which occupy a cosine band 

Ek = - J(l -I- 26n) cos(fca) - gfisB + 2JSm . (3) 

Here, k is the wave vector, m the magnetization, and a 
the distance between neighboring spins. The parameters 
and m are determined from = ^ J^^"^ cos{ka)fkdk 

and m = -\ + ^ /J"^" fkdk, where fk = (exp(efc/fcBr) + 
1)^^ is the Fermi distribution function. The magnetic 
field plays the role of the chemical potential for the JW- 
fermions. The spin thermal conductivity is given by 

Ks = / -j^ekVkts,kdk, (4) 

TT Jo dl 

where N is the number of spins per unit volume, Vk — 
h~^dek/dk is the velocity and is^k the mean free path of 
the spin excitations. The applicability of the MFT JW- 
fermion model is demonstrated by the example of the 
specific heat Cs{T,B), which is given by the right-hand 
side of Eq. ^ omitting Vkts,k from the integrand. In 
Fig. [2] (a), the calculated values of Cs{B) are compared 
with the experimental data of Ref. [12] for T = 0.35 and 
1.53 K ^21). The good agreement between the calculated 
and measured Cs{B) is obvious. Similarly, magnetostric- 
tion and thermal expansion of (C5Hi2N)2CuBr4 are well 
described by the JW fermion model [3]. 
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FIG. 2: (color online) (a) B-dependence of the magnetic spe- 
cific heat of (C5Hi2N)2CuBr4. The symbols are experimental 
data [12] and the lines are calculated within the MFT model 
(see text), (b) Comparison of the expected spin thermal con- 
ductivity Ks^min (liuc) and the experimental total thermal con- 
ductivity k" along the ladders (symbols). 



In order to compare the measured thermal conductiv- 
ity with the expected spin thermal conductivity, informa- 
tion about the mean free path Ig is required. At low tem- 
peratures, scattering by magnetic impurities and disorder 
is dominant. The scattering by phonons should be rather 
weak in the present case, because the concentration of 
phonons is much lower than the concentration of spin ex- 
citations, as the bandwidth of spinous (~ Jy ~ 3.6 K) is 
much smaller than the phonon bandwidth (of the order of 
the Debye temperature Q^:, ~ 10^ K). In previous studies 
of the heat transport in spin-l/2-ladder and spin- 
1-chain compounds |23j the impurity-limited mean free 
path is found to be (fc, T) -independent and of the or- 
der of the mean distance between impurities in the ladder 
(chain) direction. Our analysis of the low-temperature 
magnetization data of similarly grown (C5Hi2N)2CuBr4 
crystals yields an average distance of about 0.16 /im be- 
tween "defective" rungs, i.e. rungs that are not in the 
singlet state for low fields. The corresponding calcu- 
lated minimum spin thermal conductivity Ks,min(^) is 
shown for T = 0.37 K in Fig. [2] (b) and compared 
with the measured total k'L Near the critical fields, 
where the gap closes, Ks{B) should show a character- 
istic two-step increase from a negligibly small value to 
a constant of the order 100 mW m~^ K~^. Thus, one 
would expect an increase of that exceeds the mea- 
sured total K ^ Kg + Kph by more than one order of 
magnitude. Within the experimental resolution (about 
0.1 mW m-i at T = 0.37 K) our k{B) data do 

not show any indication of such steps near either Bd 
or Bc2- The absence of the expected steps cannot be 
simply explained by a very short mean free path in the 
studied crystal, because Kg < 0.1 mW m^^ K^^ would 
correspond to an unphysically small £ < 1 A. Thus, the 
thermal conductivity in (C5Hi2N)2CuBr4 is dominated 



by the phononic contribution, n Kph- This conclusion 
is true for all temperatures we studied. 

The absence of Ks most likely results from impurity 
scattering, but due to the extremely weak interladder 
coupling this scattering is not consistent with a constant 
mean free path given by the distance between impurities. 
As suggested in Ref. [23], disorder leads to a constant 
mean free path of spin excitations only in the presence 
of weak yet non-negligible interladder coupling. In par- 
ticular, the condition 



(5) 



should be satisfied, where vj- and v} are the characteristic 
velocities of, respectively, the spin excitations perpendic- 
ular and parallel to the ladder direction, and ddcf is the 
average distance between the defects. When Eq. ([5| is 
satisfied, an energy transfer between neighboring chains 
is possible and is « ddef- Otherwise, strong backscat- 
tering by impurities combined with the low probability 
of interchain transfer may lead to a spinon localization 
and thus « 0. In our recent heat transport studies in 
the spin-chain compounds copper pyrazine dinitrate and 
NENP [231 m], the left-hand side of the inequality ^ 
was an order of magnitude larger than the right-hand side 
for the temperature range where a large 1^ has been iden- 
tified. However, in the present case of (C5Hi2N)2CuBr4 
for all investigated temperatures, the estimates show that 
vj-/a < vl /ddcf- Therefore, we conclude that the absence 
of a measurable spin contribution to the heat transport 
results from a localization of the spin excitations in finite 
ladder segments of (C5Hi2N)2CuBr4. 

Theoretical efforts treating localization effects in ID 
systems have focused on the electrical conductivity [25] . 
Only recently random-disorder induced localization in 
the spin thermal conductivity of the spin- 1/2 XXZ chain 
has been specifically addressed by numerical diagonaliza- 
tion of small-size systems up to 20 spins |26J . For finite T, 
these results suggest the possibility of spinon localization 
for the easy-plane case (5 = in Eq. ([3])). Similar calcu- 
lations, which would be applicable for our case (5=1/2 
and much longer ddef), have not been published so far. 

Now we turn to another salient feature of the measured 
k{B), namely, the two minima of Kph{B) in the vicinity of 
^middle- These minima are obviously caused by the scat- 
tering of phonons by spin excitations. In order to gain 
insight into what part of the spinon spectrum contributes 
to the scattering of phonons one may use the dominant- 
phonon method. At low temperatures, when phonons are 
mainly scattered by sample boundaries (£ph = const), 
phonons with energy ~ 3.7 /cbT provide most of the total 
thermal conductivity [27] . Thus, any additional phonon 
scattering mechanism which scatters phonons around a 
certain frequency ujr, produces the strongest reduction of 
Kph at r ~ Hu!r/3.7kB. Hence, if particular spinous (in 
k space) are responsible for the minima of K(i?iiiin, T), 
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FIG. 3: (color online) (a,b) Energy of the spin excitations 
normalized by fcaT as a function of the wavevector. The dis- 
persion curves are calculated via Eq. ([3} for r and B corre- 
sponding to the minima of k{B). The arrows indicate increas- 
ing T for ka > n/2. (c) Spin excitation energy at ka = n/2 
for B — Bmin, calculated via Eq. ([3|, as a function of tem- 
perature for (C5Hi2N)2CuBr4 (BPCB) and copper pyrazine 
dinitrate (CuPzN). Both scales are normalized by the respec- 
tive exchange couplings J (see text). 



the energy of these spinons £A;(i?min, T) should scale lin- 
early with temperature. Accordingly, we calculated the 
dispersion curves of the JW fermions for the {B,T) 
values of the lower (B = -Bmini) and upper (B — -Bmin2) 
minimum via Eq. (pi and normalized them by /cbT. The 
result is shown Fig!l3l(a,b). All calculated curves Sk/k^T 
practically meet at the commensurate point ka = 7r/2, 
suggesting that the spinons with fca ~ 7r/2 are responsi- 
ble for this additional phonon scattering. 

This type of scattering seems not to be restricted to 
(C5Hi2N)2CuBr4. Minima of k{B) have also been ob- 
served in the spin- 1/2 chain material copper pyrazine 
dinitrate [H]. Fig. [sjc) compares the energies of the 
ka = 7r/2 JW-fermions e7r/2a(^min) = | - 5A*B^min + 
2J5m\, Eq. ([s]), as a function of temperature for 
(C5Hi2N)2CuBr4 and for copper pyrazine dinitrate. 
Note, that both axes have been normalized to J of the 
respective compound and that B = B ior the spin-1/2 
chains. For both materials we observe essentially the 
same behavior in the LL state suggesting a common 
and universal origin of this additional phonon scattering, 
which is most likely associated with Umklapp scattering. 
As discussed in Ref. [551 spinon Umklapp scattering in 
the presence of weak disorder can lead to minima in the 
field dependence of the spin thermal conductivity Ks(B). 
However, the minima discussed in the present paper are 
observed in the phonon part of the heat conductivity Kpjj. 

In conclusion, our measurements of the thermal con- 
ductivity in the two- leg spin-1/2 ladder compound 
(C5Hi2N)2CuBr4 provides clear evidence for the absence 
of spin- mediated heat transport. The most likely ori- 
gin is a spinon localization in finite ladder segments. 



which is favored by the extremely weak interladder 
coupling, i.e. the high degree of one-dimensionality, in 
(C5Hi2N)2CuBr4. The strongest effect observed in our 
experiments is a suppression of the phonon heat trans- 
port due to phonon-spinon scattering, which is excep- 
tionally strong close to the commensurate filling of the 
spin excitation band. 
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